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We study superradiant scattering off Bose-Einstein condensates by solving the semiclassical 
Maxwell-Schrodinger equations describing the coupled dynamics of matter-wave and optical fields. 
Taking the spatial dependence of these fields along the condensate axis into account, we are able 
to reproduce and explain many of the characteristic features observed in the experiments of Inouye 
et al. [Science 285, 571 (1999)] and Schneble et al. [Science 300, 475 (2003)], such as the shape 
of the atomic side-mode distributions for forward and backward scattering, the spatial asymmetry 
between forward and backward side modes, and the depletion of the condensate center observed for 
forward scattering. 
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The recent observation of superradiant scattering from 
Bose-Einstein condensates (BECs) 0, |3| has significantly 
extended our knowledge about collective emission pro- 
cesses. BEC superradiance differs in several key aspects 
from the well-known "conventional" superradiance in- 
volving electronically excited atoms |3)ll|- In particular, 
the role of the excited state is now played by the combina- 
tion of a BEC and an impinging laser pulse, spontaneous 
emission is replaced by spontaneous Rayleigh scattering, 
and the ground state corresponds to atoms in well-defined 
momentum side modes 0, Q ■ 

The experimental observation of BEC superradiance 
was initially interpreted in the framework of matter-wave 
stimulation 1,]. In this picture, interference between the 
recoiling atoms and the condensate at rest leads to a 
matter-wave grating, from which laser photons are scat- 
tered. This causes the matter-wave grating to grow 
rapidly in a collective, self- amplifying process. However, 
later experiments with shorter and stronger laser pulses 
led to a more profound understanding of BEC super- 
radiance and showed that the above picture was incom- 
plete. In these latter experiments, backward-scattered 
atoms were observed in addition to the forward peaks. 
As suggested by the MIT group, the new experimental 
results should not be explained in terms of optical diffrac- 
tion from a matter-wave grating, but as atomic diffrac- 
tion from the optical grating formed by the superposition 
of the impinging laser beam and the Rayleigh-scattered 
photons. An essential cornerstone of their argumentation 
was the fact that forward- and backward-scattered atoms 
show a spatial asymmetry. It was conjectured that this 
asymmetry arises because atomic diffraction is strongest 
at the edges of the condensate, where the intensity of the 
optical grating is largest. On the other hand, for pure for- 
ward scattering the picture of matter- wave stimulation is 
equally valid, which is evidenced by condensate depletion 
occurring mainly at the center where the matter-wave 
grating is supposed to be most pronounced. 

Subsequently, however, this interpretation of asymmet- 
ric scattering was questioned by Meystre and co-workers. 
According to 0, the observed asymmetry may also be 



due to the fact that larger angles are favored in back- 
ward scattering, because of a reduced energy mismatch. 
Nevertheless, their theoretical analysis was not able to 
distinguish between the two alternatives. It is one of the 
main purposes of the present paper to resolve this contro- 
versy. As our results show, the explanation put forward 
by the MIT group is indeed the correct one, thus con- 
firming the picture of optical stimulation. 

Our work is based on the semiclassical solution of the 
spatially dependent Maxwell-Schrodinger equations for 
the coupled optical and matter-wave fields. In this way, 
we are able to overcome two main limitations of the the- 
ory of Ref. 1^, i.e., the restriction to short times and 
first-order scattering (undepleted-pump approximation) 
as well as the neglect of propagation effects. So far, the 
latter have also been disregarded in most other theoret- 
ical treatments of BEC superradiance, e.g., 0, Ref- 
erences 1^ discuss spatial effects in the BEC-light in- 
teraction, but do not provide a detailed comparison to 
the results of Refs. P, Q ■ The present results show that 
the inclusion of propagation effects is essential for a com- 
prehensive theoretical understanding of BEC superradi- 
ance. Indeed, our model enables us to also reproduce 
and explain many other characteristic features, such as 
the typical momentum distribution patterns in forward 
and backward scattering and the depletion of the con- 
densate center in forward scattering. All these results 
clearly show that our model indeed captures the essen- 
tial aspects of the physics of BEC superradiance in both 
the weak- and strong-pulse regimes. Thus, whereas the 
MIT interpretation stresses the difference of the physi- 
cal pictures used to explain the observations in the two 
regimes, our approach emphasizes the existence of a uni- 
fying theoretical framework. 

In our theoretical treatment, we consider a cigar- 
shaped condensate oriented along the z axis. The 
BEC is exposed to a linearly polarized laser pulse 
£oit)ey{e'^''"=-'^'*'> +c.c.)/2, lji = ch, travelling in the 
X direction. The laser is far off-resonant from the 
atomic transition to the excited electronic state |e). 
After adiabatically eliminating the state |e), the cou- 
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pled Maxwell-Schrodinger equations of motion for the 
mean- field macroscopic wave function ipi^^ t) ^^id the 
positive- and negative-frequency components E'^^)(x, t) 
of the classical electric field read ^ [l^ 
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with S the detuning of the electronic transition, d the 
atomic dipole moment and M the atomic mass. The 
polarisation is given by P(+)(x, t) = —d\ip{x,t)\'^d ■ 
E(+)(x,i)/?i(5, P(-) = p(+)*. Note that in Eq. Q we 
neglect the external trapping potential and the atomic 
interactions, since they do not play a significant role on 
the time scales of the process under consideration. How- 
ever, it would be straightforward to include them in the 
model. 

We will solve Eqs. and (O under the slowly- varying- 
envelope approximation (SVEA). To this end, we decom- 
pose the fields as 
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E(-)(x,i) = E(+)*(x,t), with u = kc and A the av- 
erage condensate cross section perpendicular to the z 
axis. The electric field in Eq. Q contains the imping- 
ing laser pulse together with the two optical endfire 
modes, that are produced by collective Rayleigh scat- 
tering. The endfire modes, whose envelope functions are 
denoted £± {z,t), travel up and down the condensate axis, 
respectively. For simplicity, we model the applied laser 
pulse as rectangular lasting from t = Ouptoi = tf. 
In Eq. (PI, the summation is over all momentum side 
modes (n, m), with m + n even. In the side mode (n, m), 
atoms possess momentum h{nkiex + mkez) and have a 
slowly varying spatial envelope ipnm{z,t), while their ki- 
netic energy is given by hu}n,m = h'^{n'^kf + w?k'^)/2M. 
In this notation, the "side mode" (0, 0) describes the 
condensate at rest. The wave vector k is fixed by en- 
ergy conservation for the transitions between the side 
modes (0,0) and (1,±1) which initiate the process, i.e., 
Ticki — Tick ~\- huji i. Since ki — k <^ k,ki, we can ap- 
proximate U!n,m ~ (n^ + m?)ujr with the recoil frequency 
LOr = hkf /2M. We also introduce Aw = to — cui = —2ujr 
and = ujn,m-uJn±i,m±i ~ -2{±n±m + l)ujr- Note 
that in the ansatz ©-lEI we disregard the dependence 
of the envelope functions ipnm and £± on the transverse 
directions x and y. For the matter waves, this is certainly 
a good approximation since the radial degrees of freedom 
are tightly confined by the trap. For the optical fields, we 
can use this approximation, b_ecause the Fresnel number 
of the system is close to 1 [l[ 



Using the ansatz Q-Q and introducing the rescaled 
fields e± = £±/y/nuj/2eoA, Eq. (HJ in the SVEA reads 
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and the condensate length L. The first term on the right- 
hand side of Eq. (|3J) describes the quantum-mechanical 
dispersion of the envelope function, while the second one 
leads to a spatial translation with velocity v„i — mhk/M . 
The other terms describe the spatially dependent light- 
induced couplings of the momentum side mode (n, m) to 
other modes. In particular, through stimulated scatter- 
ing, an atom in a side mode (n, m) can absorb a laser 
photon and deposit it into one of the endfire modes. The 
accompanying recoil transfers the atom into one of the 
side modes (n -I- l,m ± 1). Alternatively, the atom may 
absorb an endfire-mode photon and emit it into the laser 
beam, thereby ending up in the side mode (n — 1, m± 1). 
This latter process is responsible for atomic backward 
scattering jj. Finally, following we disregard absorp- 
tion and emission between endfire modes. 

Neglecting retardation efl:ects, the envelope functions 
e± are given by 
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From these equations one sees how the build-up of 
the endfire-mode fields is driven by the coherences 

V'nm V'n-l-l, mil- 

The semiclassical Eqs. ©-0 describe the evolution 
of the system in the macroscopic regime of superradi- 
ance where the populations of atomic and optical modes 
are large compared to one. As discussed in Refs. 
such equations have to be solved with stochastic initial 
conditions (which we will call "seeds") that model the 
random character of the initial quantum noise starting 
up the process. In our case, the noise is due to sponta- 
neous Rayleigh scattering transferring condensate atoms 
into the side modes (1,±1) 12]. For any given initial 
condition, the corresponding solution of the semiclassical 
equations describes one possible realization of the exper- 
iment. The variations between different simulations re- 
flect the macroscopic effects of the initial quantum fluc- 
tuations. We have solved Eqs. for a variety of seed 
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FIG. 1: Strong-pulse regime, (a) Spatial distribution of the 
first-order forward (1, ±1) and backward ( — 1, ±1) atomic side 
modes, after applying a laser pulse of duration tf = 14 jis and 
strength g = 2 x 10^ s~^ to the condensate followed by a free 
propagation for a time tp = 25 ms. (b) Spatial distributions 
of the atomic side modes and the optical endfire modes (£±), 
at time tf. For the sake of illustration the BEG population 
(0, 0) has been divided by 4. 



functions. We find that for fixed external parameters, de- 
spite the differences in quantitative details, all solutions 
share characteristic features, such as the side-mode dis- 
tribution patterns. This shows that these features orig- 
inate from specific semiclassical dynamical mechanisms 
and not from quantum fluctuations, and it is on these 
semiclassical effects that we will focus in the following. 

In the examples shown below, we use the seed func- 
tion ipi^±i{z,0) = tpo,oiz,0)/V^ for concreteness (with 
N denoting the atom number). This corresponds to 
having one delocalized atom in each of the side modes. 
Furthermore, following the experimental data of Ref. 
0, we consider a ^'^Rb EEC with TV = 2 x 10^ and 
length L = 200 fim. The condensate is in the Thomas- 
Fermi regime, so that we can model its wave function as 
0) = ./Mz) with n{z) = C[(L/2)2 - z^]e{L/2 - 
\z\), C = 3N/AL^. The pulse duration tf and the cou- 
pling strength g are in the same regime as those used in 
the experiments (compare with Ref. 0). 

Let us first consider the regime of strong laser pulses, 
which is characterized by the superradiant gain 1] being 
much larger than the recoil frequency tOr- In Fig.^a), we 
display a snapshot of the atomic spatial distribution after 
applying a strong laser pulse to the condensate followed 
by a subsequent free propagation for a time tp ^ tf. 
Since we work with a one-dimensional model, we calcu- 
late the displacement Ax between the condensate and the 
first-order side modes in the x direction as Aa; = Vj.tp 
with the recoil velocity Vr = Tik/M = 5.9 x 10~'^m/s. 



Our result clearly reproduces the asymmetry observed in 
Fig. 2 of Ref. [21 • Moreover, we are able to verify the 
explanation for this phenomenon conjectured in As 
depicted in Fig. ^b), at time tf the atomic side modes 
and the optical field modes are indeed localized near the 
condensate edges, as suggested by the MIT group. Re- 
lated results were also obtained in a simpler model in 
Ref. 1^ disregarding backward scattering. The reason 
for the behavior shown in Fig. ^ can be inferred from the 
linearized analysis of Eqs. lEJl-lI7|) One finds that, 

during the build-up of the atomic side modes, the elec- 
tric fields £±{z,t) grow steadily in the z and —z direc- 
tion, respectively, and are strongest around the ends of 
the condensates. Equation (0 then implies that the side 
modes will predominantly grow in these areas of large 
electric fields where the atomic diffraction is largest [see 
Fig. n^b)] . The observed asymmetry is now due to the 
fact that during the subsequent free time evolution, the 
backwards scattered atoms immediately travel further 
outwards, i.e., away from the condensate center. The 
forward scattered atoms, however, will initially move in- 
wards [dashed lines in Fig.^a)]. 

Inspecting Figs. 1(A) and 2 of 01, we also find evidence 
against the interpretation of Pu, Zhang, and Meystre 
p. According to their Fig. 3(b), forward peaks should 
predominantly appear at an angle of 45° from the con- 
densate center; the experiment, however, clearly shows 
smaller angles (towards the direction of the laser pulse). 
Their Fig. 3(b) also predicts backward scattering pre- 
dominantly at 45°, but with angles larger than 45° fa- 
vored compared to those less than 45°. In the experi- 
ments there is no backward scattering at 45°, but only 
at larger angles. Furthermore, Ref. j^j predicts that the 
superradiant scattering always involves atoms with mo- 
menta q and — q simultaneously. It is thus not clear 
how atoms can appear at different angles in the forward 
and backward directions. The underlying reason for the 
failure of the model of Ref. is the neglect of spatial 
propagation effects, since our theory shows that these are 
crucial to understand the experimental observations. 

In Fig. El^a) , we show the atomic side-mode distribu- 
tion forming an X-shape pattern after applying a strong 
pulse. This figure should be compared to Fig. 1(A) of 
Ref. [3| , where the momentum distribution was obtained 
through time-of-flight imaging. In discussing this figure, 
we want to emphasize that the very appearance of an 
X-shape pattern, i.e., the suppression of off-diagonal side 
modes with |n| ^ |m|, is another independent confirma- 
tion of the fact that in the strong-pulse regime the atomic 
side modes and the optical endfire modes are located at 
the condensate edges. The off-diagonal side mode (2,0), 
for example, is resonant with the modes (1,±1), but it 
can only be populated if the (1,±1) modes overlap with 
the endfire modes £±, respectively [see Eq. (jS))]. There- 
fore, the absence of a peak for (2, 0) implies that there is 
no such overlap [see Fig. ^b)]; i.e., the modes must be 
located near the edges. The side modes (2, ±2), however, 
are readily populated since the side modes (1, ±1) over- 
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FIG. 2: Atomic side-mode distributions. Each square repre- 
sents an integrated probability pnm = J dz\ipnm{z,t)f . (a) 
Strong-pulse regime: tf = 10.6 /xs and jr — 2.6 x 10^ s~^. (b) 
Weak-pulse regime: tf = 232 /.is and jr = 5.0 x 10^ s^^. (c) 
Weak-pulse regime: tf = 291 /xs and g = 6.5 x 10^ s~^. (d) 
Spatial distribution of the condensate along the axis z corre- 
sponding to (c). 



lap with £zp . We want to stress that the X-shape pattern 
cannot be explained by models neglecting propagation 
effects, e.g., Refs. 0, Q. In these models off-diagonal 
side modes such as (±2,0) and (0, ±2) become strongly 
populated as well, so that an X pattern does not emerge. 
This again proyes the significance of spatial effects in the 
strong-pulse regime. 

We now turn to the weak-pulse regime in which only 
forward-scattered atoms with a characteristic fan-like 
side-mode pattern are observed (see Figs. 1(E)- (G) of 
[3, Fig. 1(B) of 0, Fig. 2(a) of 111). As evidenced in 
Figs. |2Ib,c), our model is able to reproduce these pat- 
terns. We would like to emphasize that the appearance 
of such patterns is a typical, not a coincidental, outcome 
in our simulations. Wc also find other types of distri- 



butions, such as the existence of population only on the 
forward diagonals n = \m\, n > 0. Such patterns, which 
have not yet been reported experimentally, appear as an 
intermediate regime between the full X-shape and the fan 
pattern. As depicted in Fig. |3d), our model also repro- 
duces the depletion of the condensate center in the weak- 
pulse regime (see inset of Fig. 1(B) of 0). Our model 
predicts that the onset of the center depletion is corre- 
lated with the appearance of the (2,0) side mode. This 
is due to the fact that the (2, 0) mode can only be popu- 
lated efficiently if the (1, ±1) modes overlap with the £± 
endfire modes, respectively. This first happens around 
the center of the atomic sample, since both atomic and 
optical modes gradually expand from the edges towards 
to the center. At that point, the condensate, from which 
the (1, ±1) modes grow, is depleted there, whereas it has 
been replenished at the edges due to Rabi flopping [T^ . 

The presented examples show that the atomic side- 
mode distributions form a sensitive probe providing de- 
tailed insight into the coupled dynamics of optical and 
matter-wave fields. A systematic investigation of these 
distributions, e.g., their dependence on laser pulse du- 
ration and strength, would constitute a form of time- 
resolved spectroscopy of BEG superradiance. 

Besides studying characteristic semiclassical effects, 
our theory also allows us to examine macroscopic quan- 
tum fluctuations in BEG superradiance. Following Refs. 
0, , we have investigated the delay fluctuations of the 
first maximum in the intensity of the emitted superradi- 
ant light. In both the weak- and strong-pulse regime the 
delay distributions are found to be well described by the 
Gumbel distribution familiar from "conventional" super- 
radiance 01 . 

In summary, we have presented an analysis of su- 
perradiant scattering from Bose-Einstein condensates in 
terms of the spatially dependent semiclassical Maxwell- 
Schrodinger equations. We have been able to reproduce 
and explain several characteristic features of the experi- 
ments, that have not been theoretically accounted for so 
far. We finally note that our theory equally well applies 
to the recently demonstrated matter-wave amplification 
|l5j | , and we expect it to be an important tool for a thor- 
ough understanding of this process as well. 
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